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Abstract 

Semi-crystalline polymers have been shown to have greatly increased thermal conductivity com¬ 
pared to amorphous bulk polymers due to effective heat conduction along the covalent bonds of 
the backbone. However, the mechanisms governing the intrinsic thermal conductivity of polymers 
remain largely unexplored as thermal transport has been studied in relatively few polymers. Here, 
we use molecular dynamics simulations to study heat transport in polynorbornene, a polymer that 
can be synthesized in semi-crystalline form using solution processing. We find that even perfectly 
crystalline polynorbornene has an exceptionally low thermal conductivity near the amorphous limit 
due to extremely strong anharmonic scattering. Our calculations show that this scattering is suf¬ 
ficiently strong to prevent the formation of propagating phonons, with heat being instead carried 
by non-propagating, delocalized vibrational modes known as diffusons. Our results demonstrate a 
mechanism for achieving intrinsically low thermal conductivity even in crystalline polymers that 
may be useful for organic thermoelectrics. 
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Bulk polymers are generally considered heat insulators due to ineffective heat transport 
across the weak van der Waals bonds linking polymer chains. However, both computational [H- 
[S] and experimental [21 IBl-ITU] studies have demonstrated that some semi-crystalline polymers 
can have large thermal conductivity exceeding that of many metals. This enhancement 
occurs when the polymer chains are highly aligned, allowing heat to preferentially transport 
along the strong covalent backbone bonds. Thermally conductive polymers could find great 
use in a variety of heat dissipation applications including electronics packaging and LEDs. 

Computational studies of heat transport in polymers have predicted the high thermal 
conductivity of crystalline polymers and also identihed key molecular features that contribute 
to this high thermal conductivity. While most studies have focused on polyethylene (PE), 
comparisons of thermal conductivity among polymers have helped to identify factors of 
particular importance in setting a polymer’s intrinsic upper limit to thermal conductivity. 
First, backbone bond strength[T] has been identified as heavily influencing group velocity, 
leading to higher thermal conductivity. Additionally, chain segment disorder [T[ [TTl |T2], or 
the random rotations of segments in a chain, has been shown to lead to lower thermal 
conductivity. 

An important goal in exploiting the high intrinsic thermal conductivity of certain poly¬ 
mers is to fabricate semi-crystalline polymers at a large scale. A morphology potentially 
suited to this purpose is the polymer brush, or an array of polymer chains attached at one 
or both ends to a substrate [13]. A promising synthesis technique, known as surface-initiated 
ring-opening metathesis polymerization (SI-ROMP), is able to uniformly grow tethered poly¬ 
mer chains from a substrate in the desired aligned structure. Polynorbornene (PNb) and 
its derivatives are well studied in the ROMP synthesis technique and are thus of interest as 
thermal interface materials [T3|. However, the intrinsic thermal transport properties of PNb 
remain unknown. 

In this Letter, we use molecular dynamics (MD) to study heat conduction in PNb. While 
PNb meets the standard criterion for high thermal conductivity, our simulations indicate 
that even perfectly crystalline PNb has an exceptionally low thermal conductivity nearly at 
the amorphous limit. We show that this low thermal conductivity arises from signihcant 
anharmonicity in PNb, causing high scattering rates that prevent the formation of phonons, 
resulting in heat being carried by non-propagating vibrations known as diffusons [T5l ITH] . 
Our work shows how intrinsically low thermal conductivity can be realized in fully crystalline 
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polymers which may be of use for organic thermoelectrics. 

We calculated the thermal conductivity of single chain and crystalline PNb using equilib¬ 
rium MD simulations with the Large-scale Atomic/Molecular Massively Parallel Simulator 
(LAMMPS). The unit cell of PNb[I^ is shown in Fig. and consists of 5-membered singly- 
bonded carbon rings each connected by a doubly-bonded carbon bridge. The simulation 
process begins by relaxing the PNb crystal in an NPT ensemble. The relaxed chain length 
obtained from the crystal is used for the isolated single chains of the same temperature. 
We use periodic boundary conditions in all directions. All simulations utilize the polymer 
consistent force held (PCFF). [T5] . 

Once the structure is thermalized, we calculate the thermal conductivity using the Green- 
Kubo formalism, 

V r°° 

= {JMJMdt ( 1 ) 

in an NVE ensemble. All simulations are run with 1 fs timesteps, and NVE ensembles are run 
for 3 ns to collect sufficient statistics for the thermal conductivity calculation. All thermal 
conductivity calculations in this work rehect the average of 5 or more identical simulations 
with different random initial conditions. Error bars rehect a single standard deviation of 
values across these simulations. 

The calculated thermal conductivity versus chain length at 300 K for both single chains 
and crystals is shown in Fig. [^. We hnd the thermal conductivity of a simple bulk PNb 
crystal at 300 K along the chain direction is 0.72 ± .05 W/mK and independent of length. 
All subsequent simulations use a chain length of 50 nm. For reference, the typical thermal 
conductivities of amorphous polymers are 0.1-0.4 W/mK[T^. The thermal conductivity 
perpendicular to the chain direction gives a thermal conductivity of 0.25 ± .03 W/mK, 
yielding an estimate of the lower limit to thermal conductivity mediated by weak inter¬ 
chain interactions. In contrast, polyethylene (PE) was found to have a thermal conductivity 
of 76.0 ± 7.2 W/mK using the same procedure, in agreement with prior works. [3] Thus, the 
thermal conductivity of crystalline PNb is very small and comparable to the amorphous 
value, in marked contrast to the large thermal conductivity of crystalline PE. 

We also calculate the temperature dependent thermal conductivity of PNb, given in 
Fig. p. This result shows a weak, positive dependence of thermal conductivity on tem¬ 
perature, a feature characteristic of amorphous materials. This temperature dependence 
contrasts with the trend for crystals, exhibited by PE, where thermal conductivity decreases 
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FIG. 1. a) Illustration of a PNb crystal, a single isolated chain, and the unit cell, consisting of 
carbon (black) and hydrogen (red) atoms, b) Thermal conductivity as a function of chain length for 
both single chains (blue triangles) and bulk crystals (red circles), indicating thermal conductivity 
is independent of length, c) Thermal conductivity as a function of temperature for PNb single 
chains and crystals. Inset: Simulated diffraction pattern of PNb crystal at 300 K along the chain 
direction with peaks demonstrating crystallinity. The temperature dependence of the thermal 
conductivity is characteristic of amorphous materials despite the crystallinity of the polymer, d) 
Thermal conductivity as a function of a dihedral energy multiplication factor for both single chain 
and bulk crystal PNb. A larger factor represents a stiffer dihedral angle. The dotted line indicates 
a factor of 1, or no change to the potential. Inset: Identical calculation for a PE crystal. PNb 
demonstrates less dependence on dihedral energy compared to PE. 
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with temperature due to phonon-phonon scattering. Although the thermal conductivity of 
PNb has an amorphous trend, it remains fully crystalline. We verify the crystallinity of the 
structure by performing simulated diffraction on several snapshots of the crystal configura¬ 
tion at 300 K. As shown in the inset of Fig. [^, clear peaks are observed, indicating that the 
structure remains crystalline even at room temperature. 

To identify the origin of PNb’s low thermal conductivity, we first consider the molecular 
features previously identified in the literature as influencing thermal conductivity. The 
first feature is, intuitively, strong backbone bonds. PNb has an entirely carbon backbone 
consisting of rings with double bonds between them. Compared to PE’s purely single-bonded 
carbon backbone, PNb is composed of equal and stronger bonds, which would suggest a large 
group velocity and thus high thermal conductivity, contrary to our calculations. 

The second criterion for high thermal conductivity is the absence of chain disorder, which 
breaks the translational symmetry of the crystal and causes scattering. PE, for example, 
has been shown to experience a dramatic drop in thermal conductivity at high temperatures 
where the repeating units begin to chaotically rotate with respect to one another. [T2] The 
onset of such disorder is connected with the rotational stiffness of the chain, which arises 
directly from the magnitude of the dihedral angle energy terms in the interatomic poten¬ 
tial. [H] Stiffer dihedral terms prevent rotation around a given bond in the molecule and 
thus inhibit segmental rotation. 

To quantitatively relate the effect of this stiffness with heat transport, we artificially alter 
the dihedral terms in the potential and observe its effect on thermal conductivity. Figure 
shows this result for both an isolated PNb chain and a PNb crystal. The inset shows the 
results of identical calculations done on PE crystals. For PNb, both single chain and bulk 
crystal show a weak dependence on stiffness. In contrast, PE shows a large increase in 
thermal conductivity as the dihedral energy is increased. These observations suggest that 
chain disorder does not play a significant role in PNb, and certainly cannot explain its 
exceptionally low thermal conductivity. This conclusion is further supported by previous 
studies of amorphous PNb that find it has highly restricted rotational movement along the 
chain. [20] 

To identify the origin of PNb’s low thermal conductivity, we next find the phonon disper¬ 
sion for PNb by calculating the spectral energy density as a function of frequency, cu, and 
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FIG. 2. a) Phonon dispersion for crystalline PNb at 300 K calculated from MD, showing contri¬ 
butions from motion in the x(blue), y(green), and z(red) directions with the chains aligned in the 
z direction, b) Zoomed in phonon dispersion from (a). Longitudinal acoustic (red) and transverse 
acoustic (blue/green) modes are clearly visible at low frequencies, c) Phonon dispersion for crys¬ 
talline polyethylene at 300 K computed identically as in (a). To avoid the dispersion being washed 
out by high intensities from low frequencies, intensities are normalized by a certain percentile of 
brightness, (a) and (b) are normalized to the 96% and (c) to the 99%. The log of the intensities 
are used in (a) and (c) to better show the whole dispersion. 


wavevector, kz, mmm 

^^{uj,kz) = CY, 


Na 


a=l 


N, 




n=l 


( 2 ) 


where G is a constant, a = {x,y,z} selects the velocity component along the specihed 
axis, n indexes the nnit cell np to Nz nnit cells in the chain, and a indexes the backbone 
atoms within a nnit cell np to Na backbone atoms. Fignre shows the intensity plots 
of kz), kz), and $ 2 ( 0 ;, kz), eqnivalent to the phonon dispersion, calcnlated for a 

PNb crystal at 300 K shown in bine, green, and red, respectively. Polymer chains are aligned 
in the z direction. The most striking featnre of the dispersion is the absence of well-defined 
modes, with a blnrred, diffnse backgronnd instead of crisp phonon bands. Fignre shows a 
zoomed in view of PNb’s dispersion, enabling clear longitndinal aconstic (LA) modes in red 
and transverse aconstic modes (TA) in bine and green to be observed below approximately 
0.25 THz. 

In contrast, the dispersion for crystalline PE at 300 K calcnlated in Fig. shows crisp 
phonon bands thronghont the freqnency range, in agreement with the literatnre [231 - I25] . The 
transition from clearly discernible modes at low freqnencies to poorly-dehned modes at high 


6 













frequencies seen in Fig. I§, for PNb has been observed before in other quasi-lD structures |26] . 
though not in polymers, and is known as the loffe-Regel crossover. [27] The crossover occurs 
when phonon mean free paths (MFPs), i, are comparable to the phonon wavelength, A. If 
the thermal vibration scatters after traveling only the order of a wavelength, the vibration 
cannot be described as a propagating wave and hence dehning a wavevector is not possible. 

To conhrm whether scattering prevents the formation of phonons, we calculated the 
phonon lifetimes and MFPs from the phonon dispersion using frequency-domain normal 
mode analysis by htting the linewidth of a chosen polarization in the dispersion to a 
Lorentzian. [211 I22| For a particular polarization, j, and wavevector, k^, the Lorentzian 
htting yields the full width at half max 2Ffc^j which is related to the mode lifetime by 

= 2Ffc^j. 

The full dispersion, ^{uj,kz), contains contributions from all polarizations, but because 
the LA modes correspond solely to atomic motion in the z direction at low frequencies, as 
evidenced by the isolated red modes in Fig. [^, kz) = kz) for frequencies less 

than 0.25 THz. All subsequent httings represent results for the LA polarization only. 

The Lorentzian htting was repeated at all wavevectors spanned by the LA polarization. 
The peaks were then ht to obtain a smooth dispersion curve, oj{kz), with an average group 
velocity Vg = 3050 m/s. To achieve a better signal to noise ratio, the values of kz) at 

closely adjacent wavevectors were averaged after matching their peaks. The data centered 
at kz = 0.02 (^) and its htting are shown in Fig. |^. 

The resulting lifetimes give the MFPs, I = VgT and are plotted, normalized by wavelength, 
A, in Fig. as a function of wavevector. The hgure clearly shows that the condition for 
lohe-Regel crossover, I ~ A, becomes satished as wavevector is increased. At higher values 
of kz than are shown in the hgure, htting the linewidth of the dispersion becomes impossible 
due to the absence of a well-dehned peak. Correspondingly, looking at the full dispersion in 
Fig. 1^, it is clear that beyond a frequency of approximately 0.25 THz, the discernible LA 
modes no longer exist and become ill-dehned. 

These calculations demonstrate that the exceptionally low thermal conductivity in crys¬ 
talline PNb is the direct result of the absence of propagating vibrations to carry heat caused 
by strong scattering. Instead, heat is carried by non-propagating, delocalized thermal vi¬ 
brations known as dihusons m [To] . This explanation is consistent with the absence of 
well-dehned modes in the dispersion as well as the characteristic amorphous temperature 
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FIG. 3. a) Spectral energy density, for the longitudinal acoustic polarization shown in 

Fig. at kz = 0.02 (^) as a function of frequency. The data points represent an average across 
nearby wavevectors and the red line represents the Lorentzian fitting that yields the lifetimes shown 
in (b). b) Mean free paths, normalized by wavelength. A, for the longitudinal acoustic modes 
of crystalline polynorbornene at 300 K as a function of wavevector. Error bars represent the 95% 
confidence interval according to the fitting. The loffe-Regel crossover occurs when £ ~ A. 


dependence observed in the thermal condnctivity of PNb despite its crystalline structnre. 

Finally, we provide additional evidence for the strong anharmonicity present in PNb by 
calcnlating the Grnneisen parameter 7 , dehned as 


d\\iu: 
5111 V 


(3) 


where u is the phonon mode frequency and V is the system volume. While 7 varies for 
different phonon modes and depends on temperature, we can estimate it by calculating 
the phonon dispersion for a range of temperatures and observing the change in volume and 
phonon frequencies. We restrict this calculation to the same LA branch used to calculate the 
lifetimes and for frequencies below the loffe-Regel crossover, where the modes are dehned. 
Three dispersions were calculated at 12 temperatures between 260 K and 400 K. The PNb 
crystals were relaxed at each temperature to obtain the system volume and dispersions were 
ht to obtain u{kz,T) over the wavevector range of the LA mode. At each wavevector kz, 
data for ln[u}{kz,T)] vs ln[I/] were linearly ht to obtain an ehective 7 . The value obtained 
for PNb after averaging over these modes is -2.8. Identical calculations for PE over its LA 










































mode yield 7 = 0.044, significantly smaller in magnitnde compared to PNb. Dne to the large 
anisotropy of the polymer crystal, it may also be usefnl to consider as the relevant effective 
volnme rather than V in Eq. The average 7 values for PNb and PE are -3.9 and -0.32, 
respectively. The difference in values by more than an order of magnitude remains. Also 
note 7 becomes negative for PE due to the negative coefficient of linear thermal expansion 
in the chain direction. The negative value for 7 for PNb is likely due to the complicated 
atomic structure and the complex molecular force field, which strongly differ from typical 
crystals where bond stretching is more commonly associated with phonon softening. These 
results strongly suggest that strong anharmonicity is present in PNb and is responsible for 
the large intrinsic scattering rates. 

In summary, we have used MD simulations to show that crystalline PNb has an exception¬ 
ally low thermal conductivity near the amorphous limit due to strong anharmonic scattering 
that prevents the formation of phonons, and thus possesses intrinsically low thermal con¬ 
ductivity. While PNb is not suitable for thermal management applications, the mechanism 
for low thermal conductivity identified here could be of interest for organic thermoelectric 
applications as perfectly crystalline polymers may have good charge carrier mobility but 
retain low thermal conductivity. This work highlights the wide range of thermal transport 
properties that can be realized in polymers, which may prove useful in employing polymers 
in thermal applications. 
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